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1. Motivation, Scope and Objectives 

In this study, design requirements for a dynamically viable four-square type gear test machine 
are investigated. Variations of four-square type gear test machines have been in use for durability 
and dynamics testing of both parallel- and cross-axis gear set. The basic layout of these machines is 
illustrated in Figure 1. The test rig is formed by two gear pairs of the same reduction ratio, a test gear 
pair and a reaction gear pair, connected to each other through shafts of certain torsional flexibility to 
form an efficient, closed-loop system. A desired level of constant torque is input to the circuit 
through mechanical (a split coupling with a torque arm) or hydraulic (an hydraulic actuator) means. 
The system is then driven at any desired speed by a small DC motor. 

The main task in hand is the isolation of the test gear pair from the reaction gear pair under 
dynamic conditions. Any disturbances originated at the reaction gear mesh might potentially travel 
to the test gearbox, altering the dynamic loading conditions of the test gear mesh, and hence, 
influencing the outcome of the durability or dynamics test. Therefore, a proper design of connecting 
structures becomes a major priority. Also equally important is the issue of how close the operating 
speed of the machine is to the resonant frequencies of the gear meshes. 

This study focuses on a detailed analysis of the current NASA Glenn Research Center gear pitting 
test machine for evaluation of its resonance and vibration isolation characteristics. A number of these 
machines as the one illustrated in Figure 2 has been used over last 30 years to establish an extensive 
database regarding the influence of the gear materials, processes, surface treatments and lubricants on 
gear durability. This study is intended to guide an optimum design of next generation test machines 
for the most desirable dynamic characteristics. The scope of this study includes: 

■ a detailed vibration excitation analysis of NASA GRC gear durability test machine, 

■ development of dynamic models for evaluation of the current test machine for resonances and 

force/vibration transmission and isolation, and ■ ; 

■ determination mechanical requirements for the connecting shaft structures and inertias for 
optimal dynamic behavior within the speed and torque ranges of operation. 
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Figure 1. A four-square type gear test machine layout. 
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This study was formed by three major tasks: (1) assessment of the gear sets of the current test 
machine, (2) development of dynamic models, and (3) parametric studies and design 
recommendations. The following sections describe each task and the studies performed to achieve 
them. 

2. Assessment of the Gear Sets of the Current GRC Test Machine 

The current GRC test machine shown in Figure 2 includes a pair of unity ratio test spur gears 
meshed in an offset configuration such that the active face width is narrower than the gear face width. 
The reaction gear set is of the same unity ratio and the center distance as the test gear pair except 
different tooth counts. As a prerequisite to any detailed dynamic analysis of the current test rigs, the 
vibration excitation from both gear pairs must be quantified. Accordingly, this task involves (i) 
inspection of a number of test and reaction gears to determine whether their nominal dimensions 
characterize them accurately, and (ii) finite-element-based contact mechanics analyses of both test 
and reaction gear pairs to determine the parameters required by the dynamic models. 

2.1. Test and Reaction Gear Tooth Surface Measurements 

This task was performed by using a state-of-the-art gear coordinate measurement machine 
(CMM) available at the Center for Research (CGR) of the University of Toledo. One pair of brand- 
new test gears, a set of brand-new reaction gears and another set of used reaction gears were 
inspected. Figures 3 and 4 show these gears while being inspected on the gear CMM. Every tooth 
on each gear was inspected in both profile and lead directions. The tooth indexing errors were 
recorded as well. 

Figures 5 to 7 illustrate CMM inspections of a typical unused test gear. In Figure 5, the profile 
traces of four representative teeth are shown. Overall quality of in profile direction is good. The 
intended profile tip relief differs slightly between the right and left flanks (for the right flank, 15- 17 
pm tip relief starting at 22 deg., and for the left flank, about 10 pm starting at 25 deg. roll angle). 
The lead charts of Figure 6 show no signs of lead errors and lead crown. Finally, the indexing 


Final Report for NASA Grant NAG3-2641 


page 4 






.0050 mm PER DIUISION 
MEASURED DEUIATION (400:1) 


DATE : 19 Mar 2002 
TIME : 00: 28 


PART # : NASA 1 


OPERATOR : M. SEHA TATLIER 


PART NAME : UNCOATEO 


SERIAL # : NASA 5 



Figure 5. Profile inspections of a both flanks of four teeth on a typical test gear. 
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Figure 7, Measured indexing errors of a typical test gear. 
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errors shown in Figure 7 are within 12 (tm suggesting that once-per-revolution errors are tightly 
controlled as well. 

A number of new and used reaction gears were also inspected. Figures 8 to 10 show CMM 
inspections of a new reaction gear. In Figure 8, the profile traces of four representative teeth indicate 
that there is a slight tip relief (about 5 jam) starting at about 15 degrees roll angle. Overall quality of 
in profile direction is exceptional. The lead charts of Figure 9 also indicate very good quality. 
Finally, the indexing errors in Figure 10 show a slight sinusoidal variation (about 15 um for this 
particular gear). For the other new and used reaction gears inspected, peak-to-peak indexing errors 
varied between 5 and 30 urn. 

In summary, both test and reaction gears measured are high precision gears with very little 
profile and lead errors. The test gear teeth have no lead crown in agreement with the nominal design 
and have a linear tip relief close to the intended design value of 15 p.m starting at a roll angle of 24.5 
degrees. The reaction gears have slight profile modifications and no lead corrections. The amounts 
of one-per-revolution indexing errors are within the range of 5 to 30 |um in diameter. From these 
inspections, it was concluded that a contact mechanics analysis of the nominal gears would represent 
the real life conditions accurately, except the influence of the indexing errors of the reaction gear 
pairs. 

2.2 Contact Mechanics Analyses 

A finite element based contact mechanics analysis of both test gears (with and without face off-set) 
were carried out to determine (i) the average gear mesh stiffness values and (ii) loaded 
transmission error excitation harmonic amplitudes. Figures 11 and 12 illustrate a contact model of a 
test gear pair with and without face offset, respectively. A face offset of 3.3 mm has been used in 
durability tests to allow four tests to be performed per gear pair. A commercially available gear 
contact analysis model CAPP was used contact mechanics analyses. Figure 13 shows the contact 
mechanics model of the reaction gear set. The design parameters for both test and reaction gear sets 
are listed in Table 1. The analyses were carried out under quasi-static conditions at 16 discrete 
rotational positions per mesh cycle. A torque value of T = 68.5 N-m (50.5 ft-lb) was used in the 
analysis. 
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Figure 8. Profile inspections of both flanks of four teeth on a new reaction gear. 
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Figure 9 . Lead inspections of both flanks of four teeth on a new reaction gear. 
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Figure 10. Measured indexing errors of a new reaction gear. 


Final Report for NASA Grant NAG3-2641 


page 13 


GEAR 10 9.2 



First, in order to establish the average mesh stiffness values, gear pairs were considered to have 
no modifications. In this case, any transmission error predicted must be due to the tooth 
deformations. The ratio of static gear mesh force applied to average value of the static transmission 
error is equal to the average gear mesh stiffness value: 

T - T 

k t - — — and k r = — 

Vp t E t r pr^r 


where k is the average gear mesh stiffness value, r p is the gear pitch radius, and £ is the average 

static transmission error. Subscripts t and r denote the test and reaction gears, respectively. Figures 
14 and 15 show the static transmission error of unmodified test and reaction gear pairs as a function 
of gear rotation at T = 68.5 N-m. From these results, the average mesh stiffness values were 
determined to be 


[8.57 x 10 7 N / m for offset gear pair 
jl . 17 X 10 8 N / m for no offset 


and 


k r =6.18xl0 8 N/m. 


The slight increase in the value of k t for the case of no face offset is primarily due to an increase in 
the contact area. It is also noted that the reaction gear mesh is nearly 7 times stiffer than the test gear 
pair with face offset. 

Next, the same test gear contact analyses were performed now for the gear pairs with the 
designed profile modifications (a 15.24 jam linear tip relief starting at a roll angle of 24.5 degrees). 
Figure 16 shows the loaded transmission error of the modified test gear pair both with and 
without the face offset. The Fast Fourier Transform (FFT) of the loaded transmission errors of 
Figure 16 yields the harmonic amplitudes of the gear mesh transmission error excitations, first three 
of which are listed in Table 2. Similarly, the transmission error of the reaction gear pair shown in 
Figure 15 was used to calculate the harmonic amplitudes of the reaction gear transmission error 
excitation that are also listed in Table 2. 
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Table 1 


Design parameters of test and reaction gears of the NASA gear durability test machine. 
(All dimensions are in mm unless specified) 


Test Reaction 


Number of teeth 28 35 

Module 3.175 2.54 

Circular pitch 9.975 7 98 

Whole depth 7.62 6.35 

Addendum 3.18 2.54 

Chordal tooth thickness ref. 4.85 3.88 

Pressure angle (deg) 20 20 

Pitch diameter 88.9 88.9 

Outside diameter 95.25 93.98 

Root fillet radius 1.02 to 1.52 

Measurement over pins 96.03 to 96.30 

Pin diameter 5.49 

Backlash ref. 0.254 0.254 

Linear tip relief 0.01524 

Roll angle of tip relief (deg) 24.5 
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Figure 11. Contact mechanics model of a test gear pair with face offset. 



Figure 12. Contact mechanics model of a test gear pair without face offset. 
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Figure 14. Static transmission error of an unmodified test gear pair without offset 
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Figure 15. Static transmission error of the reaction gear pair. 
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Table 2 


Transmission error harmonic amplitudes of test and reaction gears of the 
NASA gear durability test machine (in jum ). 


Harmonic 
Index, i 

Reaction pair 

Test pair with 
with face offset 

Test gear pair 
without face offset 

1 

0.67 

2.26 

0.56 

2 

0.42 

1.11 

0.30 

3 

0.17 

0.04 

0.13 


Final Report for NASA Grant NAG3-2641 


page 19 






3. Dynamic Models 


As stated earlier, two major objectives of the project are to (i) investigate whether any transverse 
or torsional resonance(s) exists near the operating speed of interest at the given torque value, and (ii) 
determine the stiffness properties of the for the connecting shaft structures and inertias between the 
test and reaction gear sets are suitable for isolation of the reaction gearbox from the test gearbox. 
Two dynamic models will be developed to achieve these objectives. The models will also be 
employed to define solutions if there are potential resonance and/or vibration isolation problems 
associated with the current test rig. 


3.1. Torsional Dynamic Model 

A purely torsional linear dynamic model of a NASA four-square gear durability test machine 
shown in Figure 2 is developed here as shown in Figure 17. Here, each gear of base radius r and 
polar mass moment of inertia J is assumed to be rigid except the teeth in contact. Gears are also 
assumed to vibrate in torsional direction only defined by rotational alternating displacement 6 . The 
gear mesh interface of the test gears is represented by a linear gear mesh stiffness k t and a periodic 
static transmission error excitation e t (t). Similarly, k r and e r (t) represent the mean mesh stiffness 
and the transmission error excitation of the reaction gear pair. 

3.1.1. Equations of Motion 

Considering that the subscripts rl and t2 represent test gears on shafts 1 and 2 and subscripts rl 
and r2 represent the reaction gears on shafts 1 and 2, the equations of motion of the four-degree-of- 
freedom dynamic model are given as 
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J tA\ + Ol c rM/l + r tl9tl]+ r t \k t [r tl d t{ + r t2 6 t2 ~e t (t)]+ k sX (6 [X -9 rX ) = 0 
J + r t7 c t\ r t\Qt\ + r t2^il} + r t2^t \- r t\@t\ + r t2®t2 -e t (t)]+ k s2 (9 t2 ~9 r2 ) = 0 

J r2® r2 "h ^r2^r [frl^rl *"r2^r2 ]"^ 0"2^r (/Vl^rl ^~^r2^r2 ~~ (O] — (^f2 ~9 r2 ) = 0 

Jrl#rl + r rl c r[ r rl#rl + r r20r2l+ r rlk r [r r 10rl + r r 2 0 r2 ~ e r^)]~ k sl {9 tX -0 rl ) = 0 


( 1 ) 


where c ( and c r represent the gear mesh damping values that are not shown in Figure 17 for clarity 

purposes, and k sX and k s 2 are the torsional stiffness values of the shafts connecting the test gears to 

the reaction gears. The polar mass moments of inertia of these shafts are divided into two and 
lumped to the inertias of the gears at both sides. In order to obtain a matrix form of the equations of 
motion, equation (1) can be put into the form: 


J t \6 a + r} x c t 9 tX + r tX r t2 c t 9 t2 + (r,fo + k sX )9 n + r tl r t2 k t 0 , 2 - k sl 0 rl = r tl k,e, (f) 

Jtl&ti +r t\ r t2 c t^t\ + r t2 c A2 +r t\ r t2k t Qt\ + ( r ;2^ + k s2 )0 t2 - k s2 0 r2 = r [2 k,e t (t) 

J r 2 9 r 2 ~k s2 d( 2 Jr ^r2^r^r2 "h ^r\^r2 < ~l^r\ G"l^V2^r^rl (G^t k$ 2 )9 r2 ~ ^r2^r^r(0 


(2) 


J r\9r\ -k s \6 t \ + r r X r r 2 c rQ r 2 + r } x c r 9 rX + r rX r r2 k r 9 r2 + {r} x k r + k sX )9 rX = r rX k r e r {t) 


The equations of motion are then written in matrix form as 

Mq(0 + Cq(r) + Kq(t) = F(0 

where the mass, damping and stiffness matrices are define from equations (2) as 

'j tX 0 0 0 " 

0 J t2 0 0 

M = 

0 0 J r2 0 

0 0 0 7 r i_ 


(3a) 


(3b) 
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r n c t 


C = 


r t\ r t2 c t 

0 

0 


r t\ r t 2 c t 

fat 

0 

0 


0 

0 

r r 2 c r 
r r\ r r 2 c r 


0 

0 

r r\ r r 2 c r 

r r\ c t 


K = 


r n k t +Ki 

r t\ r t 2 k t 

0 

~ k s 1 


r tl r t 2 k t 0 - k sl 

r t 2 k t + k s 2 - k s 2 0 

~ k s 2 r rt k r+ k s 2 r r\ r r 2 k r 

0 r rl r r2 k f r\k t +k sl 


The displacement vector and the external forcing vector are defined as 


(3c) 


(3d) 


Ai (o' 


r a k t e t {t) 

&t 2 M 

•, F(f) = • 

r t 2 k t e t (0 

d r 2 (0 

r r 2 k r e r (^) 

9 rl (t) 

1 

r rl k r e r (t) 


(3e,f) 


The excitation function e t (t) at the mesh of the test gear pair includes only the once/mesh 

components of the static transmission error as index errors of actual gear specimens were measured 
earlier to be insignificant 

L 

e t (0 = 2 ^ sin [ ;a V + <Pr* ] (4) 

i=l 

where e ti is the i-th harmonic amplitude of the test gear pair transmission error as listed in Table 2. 
Here, the test gear mesh frequency is defined as a> t =2nQ sl Z tl /60 where £2 5l is the rotational 
speed of the first shaft in rpm and Z fl is the number of teeth on the first test gear. As evident form 
equation (4), no once/revolution component of the test gear mesh excitation is included since the 
pitch line run-out errors, index errors and spacing errors of the actual test gears were very low as 
shown in Figure 7. 
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Meanwhile, the excitation function e r (/) at the mesh of the reaction gear pair includes not only 
the once/mesh components of the static transmission error but also once/revolution components of 
each reaction gear since indexing errors of these gears were measured to be not negligible as shown 
in Figure 10. Accordingly e r (t) is defined as 


e r (t) - ^ e r i sin i(O r f + (p n - ] + E r \ sin 


i=l 


0) r 


+ Yrl 


- E r 2 sin 


co r 


■? + Yr2 


( 5 ) 


Here e ri is the i-th harmonic mesh order amplitude of the reaction gear pair as listed in Table 2, and 
the reaction gear mesh frequency is defined as co r = 2nkl s \Z r \ / 60 where Z r \ is the number of teeth 
on the first reaction gear. In addition, the once/revolution errors of both reaction gears are included 
in harmonic form having amplitudes E r y and E r 2 and initial assembly position angles y,-i and Y r2 . 
The frequency of these two excitations are co r / Z rl and co r / Z,. 2 . 


3.1.2. Natural Modes and Forced Response 

The solution of the Eigen Value problem governing equation (3 a) yields the natural frequencies 
a>i and the corresponding mode shapes O,- ( i = 1 to 4). As the system shown in Figure 17 is a linear 
one, the forced response to the excitations of equations (4) and (5) can be calculated by using the 
Modal Summation Technique [1], In order to apply this technique, first the forcing vector F(r) of 
equation (3f) is put in the form 


F(f) : 


r t 1 

| r t2 
0 

0 


sin[i(0,r + <p„ ] + 

(=1 


0 

0 

r r2 

r r\ 


k r ^ e ri sin [i(U,.r + cp n - ] 


i=l 


+ 


' 0 ' 



' 0 

0 


CO, 


0 


■k r E r \ sin 

+ Yrl 

— • 

> 

r r l 



r r 2 

r r\ . 



M. 


k r E r 2 sin 


0) r 


■^ + Yr2 


-r 2 


(6) 
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The forced response q(Y) is then defined as 


q(0 = 


r n 



0 


n 2 
0 

L 

k i X A n -e n - sin[/(0 ; f + <p ft - + r ft - ] + 

i=l 

0 

r r2 

■k' 

1 ^ 

I 

0 



r r\ 


' 0 




0 ' 

0 



co r 


0 

r r 2 

'k r AgiE r i sin 

— t+ y rl + r £1 



_ z rl 


r r2 

M. 




r r\. 


k r ^A ri e ri sin[i(O r r + <p n - + T ri ] 


(7a) 


\k r A E 2 E r2 sin 


CO,. 

37 t+ Yr2 +T E2 


-r 2 


where the dynamic compliance matrices are 


5 

^ ti ^ 9 99 

j=i(®i ~i +2 jr£ s co s co t ) 




(7b) 






“iK 2 -/ 2 <y r 2 +2;>£>,<y r ) 


(7c) 


A£1=S- 


5=1 2 


<y, 


<W 

\ 2 


v Zrl y 


+ 


2jrCa> s CO r 


) 


'rl 


A £2 ~X' 


*s*s l 


5=1 . 2 

(<v - 


co r 


\ 2 


V Z ^y 


+ 


2yV^co^co r 


) 


'r2 


(7d) 


(7e) 


In above equations, L is the number of static transmission error harmonics that must be included in 
the analysis. Typically, it is sufficient to set L = 3 for spur gears as the higher harmonic amplitudes 
0 >3 ) are negligibly small. 5 is the total number of degrees of freedom considered in the model. 
For the system of Figure 17, 5 = 4. 
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3.2. Three-dimensional Dynamic Model 


As the torsional model proposed in the previous section should be sufficient for predicting the 
torsional natural modes and for investigating the vibration transmission/isolation issues, it is not 
capable of predicting transverse vibrations such as bearing vibrations and shaft bending vibrations. 
Such motions should become important especially when the shafts connecting the test and reaction 
gear sets are made more flexible. In this section, a transverse-torsional model of the system shown in 
Figure 1 will be developed. The model will be based on a general formulation to analyze any two 
shafts connected to each other through two separate gear meshes. This general model is to have the 
following capabilities: 

■ The model should include transverse, torsional, axial and rotational (bending) motions of the gear 
pairs that can be of spur or helical type. 

■ The model should be capable of simulating any shaft geometries including variable cross- 
sections, and the hollow shafts. The shaft transverse, torsional, axial and bending motions should 
all be included in the model. 

■ Any couplings, flywheels and rigid inertias that are mounted on the shafts including the rotatory 
inertia of the DC motor can be included. 

a The model should have the ability to included any number of rolling element bearings of any 
type supporting the shafts at specified locations positioned at any location on the shafts 

This model will be obtained by expanding an already existing geared rotor dynamics model [2]. This 
model was developed to study the coupled spur gear-shaft-bearing dynamics of a two-shaft, single- 
gear pair system and later expanded to include three shaft gear reduction units [3]. A lumped- 
parameter model of gear pair was combined with a finite element model of the shafts to predict the 
natural modes and the forced response due to static transmission error excitation. This model will be 
modified to include axial and bending motions of the shafts and gears to include helical gears, and a 
second gear pair coupling will be added to obtain the model of any four-square type closed-loop 
arrangement. The model will be made general such that any system of the same layout can be 
analyzed to allow parametric design optimization studies. In the following sections, the model 
formulation will be outlined. Only the essential details will be given here as further details can be 
found in references [2, 3]. 
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3.2.1 Overall Shaft Matrices 


The general system shown in Figure 1 has two parallel-axis shafts connected to each other at two 
places through spur or helical gears. Finite element model of each shaft can be developed by using 
Euler-Bemoelli beam elements [4]. Considering a finite element £ on one of the shafts defined by 
two nodes £ and £ + 1 as shown in Figure 18, the stiffness and mass/inertia matrices of this £ -th 
shaft (rotor) element are given as a summation of bending, torsion and axial components as 


k s( — (k s ( ) bending (k torsion + fa axial ’ (8) 

— s( ) bending +axial ( m s£~) rotation ^ ( m s(^ torsion (9) 


where the individual components of and are given in Reference [3]. These shaft element 
matrices of dimension 12 can be put into the form 



(k n )< 

(k 21 )[ 

m ,( = 

(® n )( 

(m 2 ,)[ 


_(k 2l )* 

(k 22 ) ^ 


_(m 21 )< 

(m 22 ) f 


Starting with the first shaft, if the shaft is defined by m\ finite rotor elements (£ = 1,2,..., mi ), the 
stiffness and mass matrices of the entire shaft 1 can be assembled as: 


K,! = 


( k ii)i 

(k 21 )i 

0 


(k 2 2>i + ( k ii) 2 

(k 2 1 ) 2 fa 22^2 ■*" (k 11)3 


0 0 0 
0 0 0 


symmetric 


(k22)(mi-l) + ( k u)mi 

(k 2 l)mj (k 22^m\ 


( 12 ) 
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symmetric 


( m n)i 

(m 2 i)i (m 22 )[ +(m u ) 2 

0 ( m 21 ) 2 ( m 22 ) 2 +( m ,|) 3 


0 0 0 
0 0 0 


( m 22 )(/wi-l) +( m ll)mi 

( m 2l)mi ( m 22 )m| 


( 13 ) 


The above formulation is applied to the second shaft that is described by a total of m 2 finite elements 
to obtain K j2 and M 52 • The overall shaft stiffness K 5 and mass matrices are then assembled 
as 


K 


s 


'K„ 

0 

M, = 

\M sl 

0 

0 

K,2_ 

1 3 

L 0 

M i2 . 


(14a,b) 


Both K s and are square matrices of dimension S where the total number of degrees of freedom 
is S = 6(mj + m 2 + 2) . 

3.2.2 Gear System Matrices 

A three-dimensional dynamic model of helical gear pair is shown in Figure 19 will be employed here 
[3, 5], The test gear pair is chosen as the example for this formulation formed by gears 1 1 and t2. 
Both gears are assumed to have rigid blanks that are connected to each other by a linear gear mesh 
spring k, on the plane of action in the tooth normal direction determined by the helix angle f3 t . Also 
applied in the same direction connected in series to k t is a displacement excitation in the form 
motion transmission error e t {t ) . The relative positions of the gears are such that the line connecting 
the gear centers forms the positive x-axis of the coordinate frame. In this position, the plane of action 
makes an angle \j/ f with the positive y-axis as shown in Figure 19. As the plane of action changes 
direction depending on the direction of the loading, \j/ t is defined as 
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( 15 ) 


j" <p t , 7)j : counter clockwise 

Nr. ^ : clockwise 


where <t> t is the transverse pressure angle of the gear pair and T n is the torque applied on the test 
gear t\. Helix angle (3 t is defined based on the hand of the gear t\ as 


f (+)ve if gear tl has left hand teeth 
[(-)vc if gear t\ has right hand teeth 


(16) 


Both gears are allowed to translate in x and y directions in the transverse plane and in the axial z 
direction. In addition, each gear is allowed to rotate about these three axes by Q x , d y and <9, , 

respectively. Hence, with six degrees of freedom on each gear, the gear pair t has a total of 12 
degrees of freedom that defines the coupling between the two shafts holding the gears. Equations of 
motion for gear t\ are given as: 


m t\ 'yt\ + C, p t (r) + k t p, (t) cos P, cosy/, =0 (17a-f) 

m , j x t \ + c, p, (t) + k t p, ( t ) cos P,s\ ny/, =0 

m rl z'ri - c t Pt (0 “ p, (t) sin P t = 0 

I t \9 y t\ + r n c t Pt (0 + r t\k t p t (t) sin P t cos y/ t =0 

/, A,i + r t\ c t Pt (0 + r n^t Pt W sin Pt sin V t = 0 

JtAn + r t \ c t Pt (0 + r t\k t Pt (0 cos 0, =T n 

Equations of motion for gear t2: 

m t2 ya - c t Pt (0 - Pt (0 cos p t cos y/ t = 0 (1 8a-f) 

m t2*t2 ~ c t Pt (0 ” kt Pt (0 cos Pt sm Vt ~ 0 
kCtPtif) + k t p t (Osin P, =0 
l tlQytl + r t2 c, p,(t)+ r t2 k t p t (t)sin P, cosy/, =0 
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I t2 d xt2 + r t2 c t Pt (0 + r t2 k t Pt (0 sin fi t sin \f/ t = 0 
J t2& zt2 + r t2 c t Pt (0 + r t2 k,Pt(t) COS J3 t = -T t2 


In these equations, p t (t) represents the relative displacement at the test gear mesh in the direction 
normal to contact surfaces defined by 


Pt (0 ~ [(*f l ~ x t 2 ) sin Vt + (y? 1 - 372 ) cos Vt + + r t2@zt2 ] cos Pt 

+ [”^1 +Z/2 + ( r t]&xtl + r t2&xt2) sin V, +(r n 9 ytl + r l2 0 yt2 ) cos i// t ]sin J3 t -e t {t) 


(19) 


The stiffness coupling matrix and the mass matrix of the gear pair t are obtained from equations (17- 
19) as 


( k ll)r 

( k l 2 )/ " 


M, = 

m tl 

0 

_( k 2l)r 

( k 22 )r _ 



0 

m r2 _ 


( 20 , 21 ) 


The corresponding displacement and alternating force vectors are given as 


q* = 




e t (0 


where q (l - x t i z t i 9 y t i Q x ti ® zti ] ^nd * — 1,2. 


(22, 23) 


The above formulation can be repeated for the reaction gear pair formed by gears r 1 (on the 
shaft 1) and r2 (on the shaft 2) to obtain the stiffness and mass matrices of the reaction gear pair as 
well by simply replacing the subscript t by subscript r in above equations with 

\-<p r , T, 1 '.counter clockwise 

V'r=\, T ... (24) 

\<p r , T n : clockwise 
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(+)ve if gear r\ has left hand teeth 
(-)ve if gear r 1 has right hand teeth 


( 25 ) 




Accordingly, mass and stiffness matrices and the force and displacement vectors for the reaction gear 
pair are found to be: 


'( k ll)r 

(ku)/ 

, M r = 

m rl 

0 

_(k 2 l>r 

(k 22 ) r _ 


0 

m r2 _ 


f *3 rl 1 P , , _ f(f 1 ) r } , , 

[qnj 1.(^2 )r J 


The system shown in Figure 1 has a total of two shafts and four gears forming two gear pairs 
connecting certain degrees of freedom of each shaft to each other according to the formulation given 
above. The overall gear stiffness matrix and the mass matrix, both of dimension S, can be assembled 
in the form 


( k ii), 


symm. 


o - (k n ) r 

(k 2 i), ••• 0 ••• (k 22 ) f 

0 (k 21 )r ••• 0 


( k 22)r 


(28) 


M g =Diag[--- m,] ••• m rl 


m t2 


m r2 •••] 


(29) 


The masses and inertias of other non-gear components attached to a shaft must also be included in 
equation (22) at the appropriate nodes the same way the gear masses and inertias included. 
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3.2.3 Support Stiffness Matrices 


In a typical application, each shaft is supported by least two rolling element bearings of varying 
type, size and design parameters. The most general way of describing flexibility of an individual 
bearing is to define a 6x6 stiffness matrix k bi with zero torsional terms [3, 6], If there are n b 
number of bearings in the system, an overall bearing stiffness matrix of dimension S can be 
constructed by assembling the individual stiffness matrices according to the shaft node at which each 
bearing is mounted 


0 



0 


k b l 

0 

0 


symm. 



(30) 


3.2.4 Equations of Motion 

Given the mass and stiffness matrices for shafts, bearings/case and gears, the overall mass and 
stiffness matrices of the overall system are given as 

M = Mj + M^ (31) 

K = K 5 +K 6 +K g (32) 

Finally, the force vector is defined in terms of the two static transmission error excitations as 


Final Report for NASA Grant NAG3-2641 


page 35 


F(0 = 


(fl)r«r(0 

(f 2 V/(0 

(f 2 )r«r(0 


(33) 


3.2.5 Natural Modes and Forced Response 

Equations of motion of the overall system can be written in matrix form as 

MX(r) + CX(r) + KX(0 = F(f) (34 ) 

The transverse-torsional model formulation presented up to this point did not include the derivation 
of the damping matrix C. If the damping values of each component including gear meshes, bearings 
and shafts were known, one could obtain a damping matrix that is in the same form as K. However, 
these damping values are not known in most cases. Therefore, here for practical engineering 
purposes, a set of modal damping values are used to define C instead of using a damping matrix 
that is formed by actual damping values. 

The eigen value problem governing equation (34) yield the natural frequencies co and the 
corresponding modal vectors O, (mode shapes) where i = l, 2, ...,S is the modal index. The 
Sequential Jacobi Method is used here for the Eigen Value solutions as in the case of the torsional 
model. The forced response of the system is obtained by using the same approach (Modal 
Summation Technique) as the torsional model, now with the forcing vector given by equation (25) 
where e t (t) and e r {t) are defined by equations (4) and (5), respectively. 
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4. Parametric Studies and Design Recommendations 


4.1 Dynamic Behavior of the Current Test Rig 

As the first application of the dynamic models, the current NASA GRC gear durability test 
machine is simulated. The system parameters of the test rig required by the dynamic models are 
listed in Table 3. While the torsional model of Figure 17 has only four degrees of freedom (DOF), 
the 3D model shown in Figure 20 has 246 DOF. 

The dynamic forces at the test gear mesh are of particular interest here since the durability of the 
test gear is heavily dependent on the load the gear mesh experiences. The static transmission error 
excitation amplitudes predicted by the contact mechanics models e ti and e ri listed in Table 2 are 

used in this simulation. In addition, indexing error amplitudes of E r \ = E r 2 = 10 jum are used as a 
representative value of the measured errors, and a 5 percent damping value is considered. 

Torsional and 3D model predictions of are compared in Figure 21 for the case when the test gear 
set has no face offset and no reaction gear index errors E r \ = E r 2 = 0 . The vertical axis represents 

maximum alternating force amplitude of the test gear mesh F la and the horizontal axis is the shaft 

speed in rpm. The same is presented in Figure 22 when the test gears are offset in the face direction 
as described earlier. In both figures, the predictions from both models are in reasonably good 
agreement. The differences are simply due to the factors that cannot be included in the purely 
torsional model. Comparison of torsional predictions for cases with and without offset shown in 
Figure 23 indicate that F ta are much lower when there is face off-set. 

One major observation from Figure 21 to 23 is that regardless of the model used or face off-set 
condition considered, the current test condition [7] of 10,000 rpm shaft speed is in close vicinity of 
last to natural frequencies, and hence, it is close to the primary resonance peaks. This suggests that 
the current test condition involves a significant amount of dynamic gear mesh forces in addition to 
the static force transmitted. Also observed from these figures that a minor change in the operating 
speed condition in either direction or a minor structural change to the test machine to change the 
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Table 3 


System parameters of the NASA gear durability test machine 



mass 

(kg) 

inertia 

(kg-m 2 ) 

mesh stiffness 
(N/m) 

Test gear 

0.407 

3.57 xl0~ 4 

8.57 xlO 7 

Reaction gear 

with hydraulic unit 

2.41 

2.42 xl0~ 3 

6.18x10 s 

Reaction gear 

without hydraulic unit 

3.01 

2.82xl0“ 3 

6.18x10 s 
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Dynamic Gear Mesh Force Amplitude [N] 


Torsional Model 



Shaft Speed [rpm] 


Figure 21. Dynamic force amplitude of the test gear pair as a function of shaft speed when there is 
no test gear face offset. e tl = 0.56 /Mn,e tl = 0.30 m,e r{ = 0.67 /urn, e r2 = QAljum , 
and E r \ — E r 2 — 0 . 
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Dynamic Gear Mesh Force Amplitude [N 



Shaft Speed [rpm] 


Figure 22. Dynamic force amplitude of the test gear pair as a function of shaft speed when there is 
test gear face offset. e a = 2.26 /um , e t2 = 1 . 1 1 /an , e rl = 0.67 fm , e r2 = 0.42 jum , and 
^rl = E r 2 =°- 


Final Report for NASA Grant NAG3-2641 






Dynamic Gear Mesh Force Amplitude [N] 



Shaft Speed [rpm] 


Figure 23. Comparison of test gear mesh dynamic force amplitudes for cases with and without face 
offset for E r i = E r 2 = 0, e r j =0.67 /u>n and e r 2 =0A2/jm. e t \ =2.26 fum and 
e ,2 =1.11 /m for the case of face offset, and e tl =0.56 jum and e t2 = 0.30 /urn for the 
case no face offset. 


Final Report for NASA Grant NAG3-2641 


page 42 



natural frequencies would not be sufficient to move the operating speed away from the regions of 
resonance 

4.2 Vibration Transmission-Isolation Characteristics of the Current Test Rig 

The next item investigated is the issue of force and vibration transmissibility. Ideally, the 
disturbances originated at the reaction gear pair should be prevented from traveling to the test side to 
alter the test conditions, especially gear mesh load. As the only paths of the vibration transfer, the 
connecting shafts become the focus of such an investigation. These shafts on the current machine are 
rather short and have a reasonably large diameter (L currenl =112 mm and d cunent =34 mm ) 

resulting in a relatively high torsional stiffness value ( k sl = k s 2 = 9.8 x I0 4 Nml rad ). 

In order to quantify the amount of dynamic gear mesh force caused by the disturbances of the 
reaction gear mesh, a number of limiting cases are considered. In the first case shown in Figure 24, 
only the gear mesh component of the transmission error of the reaction gear pair e r ( t ) are included 
in the simulation ( e t ( t ) = 0, E r[ -E rl = 0 ) and test gear mesh dynamic load amplitude F ta that is 
solely due to e r {t) is plotted as a function of the shaft speed. Figure 24 indicates a major resonance 
peak at 8,000 rpm exited completely by e r {t) causing Fta to reach high values. This suggests that 
the current test rig, especially its shafts, is not effective in isolating the reaction gear pair. Next, in 
addition to e r (t ) , once-per-revolution errors of the reaction gears are also included in the analysis. 
A representative value for each gear is considered E r { = E r 2 =10 /um with e t (t) = 0. Figure 25 
shows F ta as a function of the shaft speed given the relative reaction gear position angles 

Y r 2 - Yri =180° . Here, the shape the curve in Figure 24 is maintained except it is shifted upward by 
more than 300 N. In this figure, the resonance peaks are still due to e r (t) while this overall shift is 
caused by the indexing errors E ri and £ r2 . 

The total F ta due to all three excitations e t (t) , e r it) and E r \ = E r2 is shown in Figure 26 for 

three different relative reaction gear position angles y r2 - y r \ = 0> 90° , 180° . This figure compared 
to Figure 22 suggests that the influence of E rl = £ r2 is eliminated when y r2 - y, i = 0 , in order 
words, when the orientation of the high points of the indexing errors are made the same in assembly. 
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Dynamic Gear Mesh Fi 



Figure 24. Dynamic force amplitude of the test gear pair with face offset due to the e r (r) given 

e t (r) = 0 and E r \ = E r 2 = 0 . e r \ = 0.67 jam and e r2 = 0A2jum . 
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Therefore, the alternating force time history F t (t) at 10,000 rpm shown in Figure 27(a) for 
Yrl ~ Yrl = 0 has a variation at each test gear mesh cycle only. The time history in Figure 27(b) for 

Yrl ~ Yrl ~ 90° contains both mesh and shaft components superimposed on each other. The once- 
per mesh variations are caused by e,(t) and e r {t) while the once-per-revolution variation comes 

from E r i = E r 2 . Finally, in Figure 27(c) for y r 2 - y r \ — 180° (indexing errors are 180 degrees out 
of phase), the once-per-revolution variations of F t {t) are the most severe causing a significant 
chance of dynamic loading from tooth to tooth of the test gears. For instance, in Figure 27(c), 
seventh tooth of the test gears has only F ta =500 N while the 21th tooth that is 180 degrees from the 
seventh tooth has a maximum value of F ta ~ 1250 N . 

To further investigate this issue of unequal loading on each tooth, a large number of test gears 
from a previous study [8] that were run to failure were inspected on the gear CMM. A 100 percent 
inspection of each tooth were done and wear amounts on each profile were quantified as illustrated in 
Figure 28 for one of the test gears. As the wear amounts are proportional to the load applied [9], any 
unequal wear amounts would confirm the results of the dynamic model. The maximum wear 
amounts on each tooth of this test gear is plotted in Figure 29. This figure shows a sinusoidal 
variation of wear amounts from tooth to tooth reaching a maximum value of more than 150 pm at 
tooth- 12 while the wear amounts are negligible at the opposite side of the gear. Three-dimensional 
wear profiles of Figure 30 for four of the teeth (the ones marked with an arrow in Figure 29) offer 
further evidence that the once-per-revolution errors of the reaction gears were indeed not isolated to 
cause undesirable loading conditions. 

4.3 Design Parameter Studies 

Two main parameters were varied to observe their influence of the dynamic behavior of the test 
rig. These parameters are the length L and diameter d of the connecting shafts. The values of these 
parameters were varied around the current values and the resultant changes in the dynamic behavior 
were recorded. Influence of L on the torsional natural frequencies C 0 [ is shown in Figure 3 1 . Here 
a >2 and co 4 are of special interest as all of the resonance peaks in forced response curves are 
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Figure 25. Dynamic force amplitude of the test gear pair with face offset due to the 
e t (f) = 0, E rX = E r 2 =10 /urn and y r2 -y r \ ~ n ■ 0-1 = 0-67 /urn and e r2 = 
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Figure 26 Influence of reaction gear phase angles on dynamic force amplitude of the test gear pair 
for the case with face offset. e n =2.26 jum, e l2 =l.U jum, e rl =0.61jum 
e r2 = 0A2/um , E rl = E r2 = 10/Um , y rl = 0, and y r2 =Q,n/2,7t . 
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OPERATOR : M. SEHA TATLIER 


DATE : 22 Mar 2002 
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Figure 28. CMM inspection of a failed test gear. 
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Figure 28. Continued. 
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Figure 28 . Continued. 
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Figure 28. Continued. 
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Figure 28. Continued. 
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Figure 28. Continued. 
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Wear depth (pm) 




Deviation from Involute (jam) 



Roll angle (deg) 
(a) 


Figure 30. 3D tooth surface inspection of the test gear of Figure 29; (a) tooth #6, (b) tooth #12, (c) 
tooth #19, and (d) tooth #26. 
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associated with these two modes. As shown in this figure, both co 2 and ca$ are reduced as L is 
increased. We note that 0 ) 4 = 4800 Hz for L/L curren[ > 2 remaining unchanged with L. 
Considering that the mesh frequency at 10,000 rpm is f, = (10,000)(28)/60 = 4,666 Hz , an increase 
in L will not completely eliminate the resonances near the operating speed while moving them to the 

right by a certain amount as shown in Figure 32 for LI L curren[ =1,3 and 5. In this 

figure Y r 2 ~Y r \ - 0 > eliminating the influence of E r j and E r2 . The conditions for Figure 33 are the 

same as Figure 32 except y r2 -y r 1 =180°. In this figure, it is clear that, the influence of E rl and 

E r2 on F ta is minimized as the value of L is increased. In fact when L / L current = 5, the F [a 

values are the same for y r2 - y r 1 = 0 and y r2 - y rl =180° (in Figures 32 and 33). This indicates 
that a torsional shaft stiffness k s value corresponding to Li L currenl =5 is sufficient to isolate the 
excitations created by the reaction gear pair from the test gear pair. Time histories F, ( t ) shown in 
Figure 34 for L / L current = 1, 3 and 5 support this conclusion further. 

Finally, the influence of shaft diameter d on the torsional natural frequencies &»,• is shown in 
Figure 35. Here ty 4 does not change significantly when d/d curren[ < 1 . Meanwhile, increasing 
beyond d curren[ does not improve on the resonance conditions as co 2 is also brought near the 
operating speed range as well. In the process, the vibration isolation conditions are also worsened. 

One other potential design change is to add inertias (flywheels) between the gear sets to further 
improve the dynamic isolation while moving the resonances away from the operation speed point. 
As the point of application of the flywheels become important, the 3D model is used here for this 
comparison. Figure 36 illustrates the 3D dynamic model of the NASA gear durability test rig with 
LI L current = 5 . In Figure 36(a), there is no added inertia while two inertias (one on each shaft) are 
included near the test gears in Figure 36(b). These flywheels are selected to have the same polar 
moments of inertia as the reaction gears, i.e. J ^ = J ri = 2.8 x 10 kg-m . In order to prevent any 

increase in the translational vibrations of the shafts, another set of bearings are mounted on the other 
side of the flywheels. In Figure 37, the resultant dynamic responses of these two configurations 
shown in Figure 36 are compared to that of the current system shown in Figure 20. Here, added shaft 
length alone is not sufficient to move the resonance peaks away from the operating speed of 10 , 000 . 
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^ current 


Figure 31 


Influence of the length of the connecting shaft L on natural frequencies 
A current — nun and d curreni — 34 171171 . 
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Dynamic Gear Mesh Force Amplitude [N] 



Shaft Speed [rpm] 


Figure 32. Influence of shaft length L on the dynamic force amplitude of the test gear pair with face 
offset. e ti ~ 2.26 jum , e t2 =1.11 jum , e r[ = 0.67 Jim , e r2 =0A2/um, E rl - E r2 - lOjim , 
Yrl “ Y r2 “ ® > L current ~~ 1 12/WW and ^ current = 34 A7Z/77 . 
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Dynamic Gear Mesh Force Ampli 



Shaft Speed [rpm] 


Figure 33. Influence of shaft length L on the dynamic force amplitude of the test gear pair with face 
offset. e tl = 2.26 jum, e t2 =l.ll jum y e rl = 0.67 jum , e r2 = 0.42 /um , 
Erl ~ E r 2 — 10/Mn , y r | —0, Y r2 —7l % E curren t 112 /72/T2 and d current = 34-t?i?n. 
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Mesh Cycles 


Figure 34. Dynamic force time history of the test gear pair with face offset at 10,000 rpm shaft 
speed. e t i=2.26pm, 2^2 =1.11 pm, e r ^ -0.61 pm e r2 =0.42, E r \ = E r 2 =10 pm , 
Y r 1 ~ 0 and Y rl — ft • (<0 L! ^current — 1 , (b) L! L currenl — 3 , and (c) L,! E currenl — 5 
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(b) L/L 



Figure 34. Continued. 
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Figure 34. Continued. 
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d / d current 


Figure 35. Influence of the diameter d of the connecting shaft on natural frequencies of torsional 

system. L current = 1 12 mm , d current =34 mm . 
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Figure 36. 


3D Model of the test rig for LI L currem = 5 ;(a) without any additional inertias, and 
(b) with additional inertias near test gears. 
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rpm. Meanwhile, the configuration shown in Figure 36(b) with added inertias eliminated most of the resonance 
peaks around 10,000 rpm. This suggests that a combination of added inertias and added compliance between 
the test and reaction gears can improve the resonance and vibration transmission characteristics of the test 
machine significantly 

5. Conclusions 


5.1 Summary 

In this study, dynamic behavior of the NASA gear pitting machine is investigated. The current 
test and reaction gears are measured using a precision gear CMM and analyzed using a gear contact 
analysis model for describing their geometric and elastic properties and excitation mechanisms. Two 
different dynamic models, a purely torsional mode! and a three-dimensional model, were developed 
to predict the dynamic loads experienced by the gear meshes. The interactions of test and reaction 
gear pair were investigated. A number of potential design options were studied to quantify their 
impact on the dynamic and force isolation behavior. 

5.2 Conclusions and Recommendations 

Based on the results of this investigation, it can be concluded that the current test machine as 
operated at the current test conditions experiences dynamic loads that are due to the combined 
contributions of the excitations from both test and reaction gear sets. A number of resonance peaks 
exist near the operating test speed condition. The connecting shafts are too rigid to isolate any of the 
excitations from the reaction gear set from the test gear set. Accordingly, for the next-generation test 
machines, following changes to the current test rig are recommended: 

■ The overall flexibility of the shafts must be increased by changing the dimensions of the 
shafts. Results of this study indicates that at the current diameter, the distance between the 
test and reaction gear pairs should be at least five times the current dimension for an 
acceptable isolation of the reaction gear pair from the test gear pair. If the shafts to be 
maintained at the current length, torsionally flexible couplings can also be used for this 
purpose. 
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■ Additional flywheels of certain amount of inertias can be implemented near the test gear 
pair to move the resonance peaks away from the operating speed condition resulting in a 
relatively resonance-free dynamic condition. 

■ An in-line torque meter should be incorporated to monitor loads experienced by the test 
gear set continuously. 

■ The face offset of the test gears should be eliminated and a slight lead crown should be 
incorporated in test gear design in order to improve the contact patterns and eliminate 
asymmetric (twisting type) tooth deflections that result in significant edge loading. 

■ The dynamic behavior of the new test rig layout must be studied using the dynamic models 
developed to prevent any undesirable conditions due to flexural motions of shafts and 
bearings. 

For the future use of the current machine, following operational changes to the current test procedure 
can be recommended: 

■ The face offset of the test gears should be eliminated for the same reason as above. 

■ The reaction gears should be paired based on their indexing errors so that each have the 
same amount, if possible. In addition, they should be assembled in an “in-phase” condition 
for a cancellation of their influence on the transmission error excitation. 
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